Local Moment Formation of an Anderson Impurity on Graphene 
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We study the property of a magnetic impurity on a single-layer graphene within an Anderson 
impurity model. Due to the vanishing local density of states at the Fermi level in graphene, the 
impurity spin cannot be effectively screened out. Treating the problem within the Gutzwiller ap- 
proximation, we found a region in the parameter space of U-E^ where the impurity is in the local 
moment state, which is characterized by a zero effective hybridization between the bath electron 
and the magnetic impurity. Here U is the onsite Coulomb repulsion of the impurity electrons and 
Ef the impurity energy level. The competition between U and E-^ is also discussed. While larger 
U reduces double occupation and favors local moment formation, a deeper impurity level prefers 
double occupation and a nonzero hybridization and thus a Kondo screened state. For a fixed U, by 
continuously lowering the impurity level, the impurity first enters from a Kondo screened state into 
a local moment state and then departs from this state and re-enters into the Kondo screened state. 

PACS numbers: 73.20.Hb, 71.27.+a, 75.30.Hx 
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I. INTRODUCTION 

The research on a single layer of graphite - graphene - 
has been advanced tremendously after its successful ex- 
foliation,— the importance of which cannot be overstated. 
Due to its exceptional electronic, thermal, and mechan- 
ical properties, graphene virtually finds its potential in 
every aspect of modern technology. For example, it could 
have tremendous applications in microelectronics due to 
its high electron mobility. It also has a great deal of 
potential in spintronicS)^'^ where the spin current is uti- 
lized. From the theoretical perspective, it has been a 
test ground for many fundamental physics laws due to 
its similarity to massless Dirac fermionSf^i^ where the 
"speed of light" is replaced by the graphene Fermi ve- 
locity Vf ^ lO^^Vs. 

The property of a magnetic impurity in a metallic host 
is an important problem in condensed matter physics. In 
a conventional metal, the Kondo effect, where the itin- 
erant electrons completely screen the magnetic moment 
of the impurity, is always present at low temperatures. 
The impurity property is not sensitive to the detailed 
band structure of the host metal. The most important 
element from the impurity bath is its density of states 
(DOS) at the Fermi level, p{u! — 0), and it enters the 
physics through the scattering rate F = p{0)V, with V 
being the hybridization matrix between the impurity and 
host electrons. This is due to the fact that in conventional 
metals, the DOS around the Fermi energy is almost a fi- 
nite constant. However, it has been theoretically found 
that vanishing host DOS (power law or otherwise) could 
lead to significant changes to the Kondo physics The 
main finding, among many others, is that the small avail- 
able states around = for conduction electrons makes 
the effective hybridization between the band electron and 
the impurity vanishingly small. This ineffective screening 
of the impurity moment causes the impurity to decouple 
from the conduction media and form a free local moment 



(LM). Withoff and Fradkiii^ found that the Kondo cou- 
pling, J, has to be larger than some critial value for the 
Kondo effect to take place. However, the dependence of 
the magnetic impurity state on the conduction electron 
DOS is a tricky one. For exmaple, numerical renormal- 
ization group calculation shows that the impurity state 
depends not only on the index of a power law DOS but 
also on the size of this power law region)^ 

The hexagonal structure of graphene has two inequiv- 
alent carbon atoms, A and B, in one unit cell. In the 
tight bind approximation, the graphene electronic struc- 
ture is modelled by a nearest neighbor hopping integral 
t. Further neighboring hopping integrals are small and 
are usually neglected. The conduction and valence bands 
then touch at the corners of the graphene Brillouin zone, 
K and K', which are called the Dirac points. The disper- 
sion near these Dirac points is linear in k and the DOS 
is also linear in energy. For all practical purposes, this 
tight binding description is quite accurate. Many pecu- 
liar properties of graphene are closely related to this lin- 
ear dispersion and the corresponding linear DOS at low 
energies The state of a magnetic impurity on graphene 
is expected to be very different from that on an ordi- 
nary metal. However, the existence of two inequiva- 
lent Dirac points could also complicate the simple lin- 
ear DOS argument since this may lead to multi-channel 
Kondo effect liSi The specific problem of a magnetic im- 
purity on graphene has recently been studied by several 
groupsjiir— Besides its theoretical importance as a pro- 
totypical many-body toy model, it is also important to 
understand magnetic impurities on graphene from the 
application perspective. It was observed that introduc- 
tion of vacancies into graphite or graphene could indeed 
induce ferromagnetismJ^"— 

In this paper, we explore the properties of a magnetic 
impurity in a single-layer graphene within an Anderson 
impurity model. Within the Gutzwiller approximation, 
we are able to obtain a phase diagram within the parame- 
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ter space spanned by the impurity level and on-site Hub- 
bard repulsion. We identify a parameter region, where 
local moment is formed on the impurity. In addition, a 
re-entrant Kondo resonance behavior is found. The effect 
of graphene chemical potential, fj., is also discussed. 

The outline of the paper is as follows: In Sec. |lll 
the model Hamiltonian for an Anderson impurity in the 
graphene is introduced. The Gutzwiller procedure to- 
ward the solution is explained. In Sec. lIIIl we present the 
results on the local moment formation and Kondo reso- 
nance state as a result of the competition between the 
impurity level and on-site Hubbard repulsion. Finally, a 
summary is given in Sec. IIVI 

II. MODEL HAMILTONIAN AND 
THEORETICAL METHOD 

We start with a single-impurity Anderson model in 
graphene. The Hamiltonian is partitioned into three 
parts — the pristine graphene subsystem Hq, the impu- 
rity subsystem i^imp, and the hybridization -ffhyb between 
these two subsystems, that is, 

v. = Hq+ iJimp + -ffhyb , (1) 

where the three parts can be written out explicitly, 

^o = -*E^('*)°kA- + H.c. , (2a) 

kcr 

CJ 

^hyb=^E«L/- + H-C-- (2c) 

We take a k-independent hybridization between car- 
bon TT-electrons and the impurity spin, Vk ~ V. Here 
^(k) = X]i=i ^''"'S where iji, i = 1,2,3, are the vectors 
connecting to the three nearest neighbors of a carbon site, 
and Af is the number of unit cells. As mentioned above, t 
is the nearest neighbor hopping integral, t = 2.7eV. In the 
rest of the paper we take t as the unit of energy. The im- 
purity is taken to be at the origin of the A-sublattice. For 
the more general case where a magnetic adatom is consid- 
ered, the impurity could be anywhere and hybridization 
between the impurity and its all neighboring sites should 
be considered. As discussed by Uchoa et al.^ this leads 
to an effective k-dependent hybridization l^(k). 

The rich physics of the Anderson impurity model 
comes from the on-site Hubbard interaction U on the 
impurity. For small U, a Hartree-Fock approximation is 
usually sufficient. For the cases of intermediate to large 
U, more degrees of the correlation effect must be taken 
into account. Here we take a projection approach and 
treat U in the Gutzwiller approximation. In this approx- 
imation, the effect of U is taken into account by reducing 
the effective hybridization between the conduction elec- 
trons and the impurity electrons. For example, in the 



U = 00 limit, the process /Jaoo- is prohibited if there is 
already a a-f electron at the impurity site, while there 
is no such restriction if [/ = 0; see Re^ for a detailed 
discussion. On the mean field level, we renormalize the 
hybridization strength V with Va- — g^V , < 1, and re- 
place the interaction term with Ud, with d = (/|/t/|/4,). 
Here g^, is the Gutzwiller factor that depends on the im- 
purity number occupation (n-^ = n| -|- being the 
total impurity occupation) and double occupation d. The 
spin dependent Gutzwiller factor is, 

Jil-nf +d){nl-d) + Jd{ni~d) 
5<T = , ■ (3) 

Note that this form of g„ explicitly breaks the spin ro- 
tational invariance, which the model Hamiltonian pos- 
sesses. However, in the following, we will only concern 
ourselves with the stability of the nonmagnetic Kondo 
screened state, = n^, and thus concentrate on the 
calculation where 5t ~ 54- ~ 5- ^ constraint term, 
X^cr ^T^'^i ~ '"^a) ^ '^^ '^^^'^ added to the mean field Hamilto- 
nian. Here the Langrange multipliers, X„, serve to renor- 
malize the impurity energy level. Although this term is 
formally the same to a direct Hartree decoupling of the 
interaction term, lining, the self-consistency condi- 
tions are quite different, and thus would lead to differ- 
ent physics. It is noted that this approach is equivalent 
to the Kotliar-Ruckenstein slave-boson approximation;^ 
which has been successfully applied in the study of An- 
derson impurity modelsi^ After the above mean field 
procedure, we have the renormalized hybridization and 
impurity Hamiltonians, 

^hyb = -^^KaL/.+H.c., (4) 

kcr 

H;n,p = J2{Ef + X,)flU + Ud. (5) 

a 

As such, the effective non-interacting Hamiltonian is, 
"H = Hq + -ffhyb + ^imp- Thcrc are two major differences 
comparred to the Hartree-Fock mean field approxima- 
tion. First, the hybridization strength V is renormalized 
with a factor g^ that is dependent on the impurity oc- 
cupation 71^ and double occupation d, which is in turn 
dependent on U. The second difference is that the im- 
purity level is renormalized by the Lagrange multipliers 
Act , instead of Un^. In the Hartree-Fock approximation, 
the impurity level renormalization U rig becomes unphys- 
ically large for large U . 

In the limit ?7 — >■ oo, the double occupation d should 
tend to zero. By looking at the form of ga in Eq. ([3]), 
it is clear that if at the same time goes to 1, then 
ga approaches zero. In this limit, the magnetic impurity 
and the conduction band decouples since the effective 
hybridization = QaV is now zero. The impurity is 
in the LM state. It is interesting to realize that this is 
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analogous to the Brinkman-Rice type of metal-insulator 
transition in the Hubbard model at half filling)^ except 
that here one requires "half filling" on the impurity site. 
Unlike the half filled Hubbard model, in addition to the 
chemical potential /i, the impurity number occupation 
is also dependent on several other aspects. It obviously 
is a function of the impurity level . It also depends 
on U, since the interaction modifies the impurity level 
through Aq.. The conduction electronic structure and the 
hybridization V also affect the impurity occupation. One 
can then tune these parameters so that the impurity site 
is half filled, and then expects to have an LM state for 
appropriate values of U, where d = 0. It is also worth 
pointing out that having a state with = 1 does not 
mean that the impurity is in the LM state. The criterion 
is to have a zero . This will become apparent when we 
present the results below. 




FIG. 1: (Color online) The Gutzwiller factor g„ (a) and the 
effective impurity level + Act (b) as a function of the on- 
site interaction U for two different chemical potential values 
= and ^ — —0.13. The impurity level is fixed at E-^ = —2 
and the hybridization strength is set at V = 0.52. 

We now set out to solve the effective Hamitonian H. 
First of all, 'H is now non-interacting and can be readily 
solved by the usual Green function method, except that 
we need to take care of the self-consistency equations 
of d and Ao-, which are obtained by minimizing the free 
energy. At zero temperature, T = (on which the rest 
of the calculation is based), the free energy equals to the 
ground state energy and the self-consistency equations 
are obtained by taking partial derivatives of (H) with 
respect to d and n^: 



^ + ^E^[(«U) + (/]«o.)]=o, 



(6) 
(7) 



where the expectation values (■ ■ ■ ) are taken in the ef- 
fective Hamiltonian T-l and are determined from their re- 
spective retarded Green functions. Applying the equa- 
tion of motion method, we find the following retarded 



Green functions. 



-Ef -X^-V^^G^iu)- 



(8a) 



(^+)2-|e(k)|2 UJ. 



KG^ ,t(w), (8b) 



where uj~^ = uj + iO~^ and G^{uj) is the retarded Green 
function for the clean graphene: 



(9) 



and 



The expectation values in Eqs. ^ and ([T]) 

= {fafij) are then computed as (O1O2) = 
-i/dw luiG^^Q^{uj)fFD[<^ - /i), with fpD the Fermi- 
Dirac distribution function. For a given set of parame- 
ters, {U,^E^ , V, ^}, Eqs. jH) and ^ together with the 
form of H yield a set of nonlinear equations of d, and 
Act, which can be readily solved. The impurity DOS is 
given by pl{u) = -ifm .{uj). 

J a- -ij a- 



III. RESULTS AND DISCUSSIONS 

Let us first fix the impurity level E-f and see the ef- 
fect of U on the impurity state. Figure [IJa) shows the 
Gutzwiller factor as a function of U for two different 
cases: fi = and fi = —0.13. The impurity level is at 
E^ = — 2. For small U the Gutzwiller factor g„ is unity, 
which is expected, since in the non-interacting impurity 
situation, there should be no renormalization of the hy- 
bridization strength. As U increases, g^r continuously 
decreases. There is a critical value oi U, Uc, in each case 
where go- is zero. At this point, the effective hybridization 
between the conduction electron and magnetic impurity 
is zero and an LM state forms. As discussed earlier, the 
impurity occupation has to be = 1/2 in the LM state 
(see Fig. [2I^a)). The formation of the impurity LM state 
can be envisioned from two facts: g^ approaches zero 
and the effective impurity level {E-^ + Ac) — > fi as U 
goes toward Uc (see Fig. [IJb)). From the Green function 
Eq. ([5a|) . it is then clear that, if Va ^ at the same 
time, the impurity DOS is narrowed to a level just at the 
Fermi energy. That is, the decoupling from the conduc- 
tion band and the effective impurity level reaching the 
Fermi energy have to happen at the same time for an 
LM state to be realized. The LM state is stable against 
further increase of U since the impurity is already decou- 
pled from the conduction band. Larger U will not change 
the double occupation, which is already at its lowest pos- 
sible value d = 0, and thus will not change the impurity 
occupation as can be seen from Eq. ([7]). Note that before 
the system reaches the LM state, the impurity occupa- 
tion does depend on U by shifting the effective impurity 
level through Ao- and in this regime d is still dependent 
on U. 
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FIG. 2: (Color online) (a) The impurity occupation as a func- 
tion of the onsite interaction U for ^ = and /i = —0.13. In 
both cases, the spin-resolved occupation number goes to 1/2 
as U approaches the critical value, (b) The resonance peak 
around /x for several values of U. The resonance peak moves 
progressively toward fi and the peak becomes sharper and 
more symmetric as U is increased. The impurity level is fixed 
at E-^ = —2 and the hybridization strength is set atV = 0.52. 



The different critical value Uc for the two cases, /i = 
and yU = —0.13 suggests that graphene DOS at the Fermi 
level plays a central role. The graphene DOS is zero at 
/J, = 0, and the hybridization between the host and im- 
purity is inefficient to fully screen out the impurity spin. 
By comparison, for the fi = —0.13 case, where there is 
a finite DOS at the Fermi level, a larger value of Uc is 
needed to bring the impurity into the LM state. We can 
look at how the resonance peak evolve as we ramp up U. 
As shown in Fig. [^IJb), the resonance peak moves toward 
the Fermi level and also gets sharper for increasing U. 
Due to the singularity of the graphene DOS at w = ±<, 
the impurity DOS also has peaks around these energies 
(not shown in figure). But as U increases, the weight 
of these states diminishes and the impurity occupation 
mainly comes from around the resonance energy. As the 
resonance peaks move from below toward the Fermi en- 
ergy, they also become more symmetric, suggesting the 
effect of V^Goo{uj) term in Eq. ((8a|) falls down. 

Let us now look at the effect of the impurity level. Un- 
like the case of varying U, the LM state is not stable if 
one lowers the impurity level past a critical value. Fig- 
ure[3Ua) shows the impurity occupation as a function of 
the impurity level for three different values of U. We 
can see that as the impurity level goes deeper inside the 
conduction Fermi sea, the occupation increases. It is in- 
teresting to notice that there is a fiat region of the curve, 
where the impurity occupation seems independent of E-^ 
for moderate on-site Hubbard interaction U, for example 
U = 2.6 in Fig. ^a,). The impurity occupation in this 
region is 1, which suggests that an LM state has possi- 
bly formed. For smaller U, the flat region shrinks and 
is slightly dependent on E-f , for example, U = 2.3 in 
Fig. EJa). However, although the curves in Fig.[3I^a) look 
similar, the physics are fundamentally different for these 




FIG. 3: (Color online) (a) Impurity occupation as a function 
of E-^ for three different values of U. (b) The Gutzwiller 
factor (7o- as a function of E-^ for the same three C/'s as in (a); 
V = 0.52, fi = -0.13 in both plots. 



different values of U. The large U value curve represents 
a transition from a Kondo screened state to an LM state 
and then back to a Kondo state as the impurity level 
is further lowered, while there is no such transition for 
the smaller U curve. This can be more convincingly seen 
from Fig.[2Jb) where is plotted against E^ for the same 
values of U as in Fig. [3{a). For U = 2.6, the Gutzwiller 
factor decreases as E-^ goes down, reaching zero at 
the upper critical value E^^ = —1.1. Before E^ reaches 

the lower critical value iS/j — —1.76, the only solution to 
the system is the decoupled conduction band and impu- 
rity. The region between E^^ and E^2 ^^so corresponds 
to the flat region where the impurity occupation is 1 in 
Fig. [3Ja). The upper critical value Ej^^ and lower critical 
value E^^ of Ef are identified from the zero segment in 
ga , which is much clearer than from the fiat region of im- 
purity occupation as in Fig. ^a.) . One can also see from 
Fig. El^b) that smaller U does not produce an LM state 
for all valules of E-^ : The Gutzwiller factor is always fi- 
nite for U = 2.3. There is a critical value of U such that 
below this value there is no LM state for any value of 
Ef. This is identified as when e[^ = E^^ = -1.35 for 
U = 2.44. Comparing Fig. [S^a) and Fig. ^h), we see 
that g^r ~ is BL more natural and physically clear way of 
identifying the LM state. 

The entrance and departure of the LM state upon low- 
ering the impurity level suggests the competition between 
U and E^ . Large U would inevitably put the impurity in 
the no-double-occupancy state, d = 0. In Fig. 2] we plot 
the double occupation, d, against E-^ . When E^ is close 
to zero, due to small available graphene electron DOS, 
the impurity occupation is small and thus the double oc- 
cupation is also small, since d < n^. Lowering E-^ will 
increase d by increasing the overall impurity occupation 

. But eventually the effect of U kicks in and further 
lowering E^ actually makes double occupation unfavor- 
able. As approaches 1 and d goes to 0, g^r approaches 
0; where we have an LM state. But further lowering E-f 
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FIG. 4: (Color online) Impurity double occupation d as a 
function of the impurity level for the same three values of U 
as in Fig.O Inset: Zoom in of the main figure for from -2 
to 0. 



favors a state, where is larger than 1. This inevitably 
allows double occupation and moves the impurity out of 
the LM state and back to the Kondo screened state again. 
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FIG. 5: (Color online) The regions of LM state and Kondo 
screened state on the U-E^ plane for different sets of V and 

Finally, we put together the effects of U and on 
the impurity state and map out the regions of Kondo 
screened state and LM state on the U-E^ plane. We 
plot the "phase diagram" of the impurity state in Fig. [5] 
for three sets of {ji^V). The general feature is that for 
moderate J7, as we lower the impurity level, the impurity 
goes from the Kondo screened state to a decoupled LM 
state, and stays in this state until the level is low enough. 
At this point, the impurity and the host graphene start to 
talk to each other again. However, if we look along the U - 
axis, after the impurity enters the LM state, it remains in 
that state even if we further increase U . Comparing the 
sets of curves, we see that larger hybridization strength 
favors the Kondo screened state and thus gives a smaller 
LM region. Also, moving away from the Dirac points 



makes it less susceptible to LM formation. 



IV. SUMMARY AND CONCLUDING 
REMARKS 

In conclusion, based on the Gutz wilier approximation, 
we have carried out a detailed study on the LM forma- 
tion of a single magnetic impurity in graphene. Due to 
its linear DOS near the Fermi level, graphene allows the 
possibility of LM formation. It has been shown that for 
moderate impurity level and large on-site Hubbard in- 
teraction on the impurity, an LM state is favored over a 
Kondo singlet state. We have explored the competition 
of impurity level and on-site Hubbard interaction and 
mapped out the LM region on the U-E^ plane. We em- 
phasize the importance of identifying the LM state by the 
vanishing of the effective hybridization as opposed to the 
unit impurity occupation, n-f = 1. The effect of graphene 
Fermi level is also discussed. It is remarkable that the 
simple Gutzwiller approximation captures the essential 
physics in this model. The Gutzwiller approximation be- 
ing a mean field approximation of the projection method 
does not handle the excited states very well. But as we 
have seen that the spetral weight of the impurity state 
is predominantly centered around the resonance energy. 
This feature justifies the applicability of the method used 
in this study. 

The results presented here for a single Anderson im- 
purity also remind us of the so called Kondo breakdown 
(KB) phenomen a^^'^^ in a lattice model, where a heavy / 
band decouples from a wider c band. In the KB regime, 
the decoupled / band would usually be polarized due to 
exchange interaction among the /-electrons. The mag- 
netic state of the LM impurity is not determined in the 
single impurity case. A small magnetic field will align it 
in the direction of the field, due to the Zeeman energy. 
For more than two impurities, RKKY interaction occurs, 
and the magnetic state of the impurities in the LM state 
depends on the nature of this interaction. For example, 
for two impurities situated at the lattice sites the RKKY 
interaction could be ferromagnetic or antiferromagnetic, 
depending on whether they are on the same or opposite 
sublattice.i^"^ The sign and strength of the RKKY in- 
teraction also depend on the graphene doping level, and 
whether the impurities are away from the carbon atoms. 
In such circumstances, the real-space Gutzwiller approx- 
imation is an ideal tool for investigating the impurity 
magnetic states. 
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